Computer  Science  Department 


TECHNICAL  REPORT 


FINDING  EULER  TOURS  IN  PARALLEL 
By 
Mikhail  Atallah 
and 
Uzi  Vishkin 

Technical  Report  #134 
Ultranote  #73 


(a) 


o 


n 
I 

Eh 


03 
u 
D 
O 
-P 

u 
0) 


U    3  . 

W  -C  u  M 

O  -H  C  M 

\4->  ^a  vj 

:d  <c  c  m 


NEW  YORK  UNIVERSITY 


Department  of  Computer  Science 
Courant  Institute  of  Mathematical  Sciences 

251  MERCER  STREET,  NEW  YORK,  N.Y.  10012 


W"^^         «S      -  ^         ^*         M         ^  M         •  • 

:#  ':#  ^in  0  -^li' ' «   ii  .<•    ii    ii    4 
'  «^  -Hi'   ip'   «   li    #    i»    «^    »    «i 
-fi  -0   0 "  0   IB    w    ii    <ii    w    4 


FINDING  EULER  TOURS  IN  PARALLEL 
By 
Mikhail  Atallah 
and 
Uzi  Vishkin 

Technical  Report  #134 
Ultranote  #73 


lep^l 


Finding  Euler  Tours  in  Parallel^'*' 
November  1982.  First  revision:  February  1983.  Second  revision:  May  19RA, 


Mikhail  Atallah 
Department  of  Computer  Science 
Purdue  University 
West  Lafayette,  Indiana  47907 

and 

Uzl  Vlshkln^^) 

Department  of  Computer  Science 

Courant  Institute  of  Mathematical  Sciences 

New  York  University 

251  Mercer  St.,  New  York,  NY  10012 

Abstract. 

The  problem  of  finding  Euler  tours  in  directed  and  undirected  Euler  graphs  Is 
considered.  We  give  O(log  |v|)  time  algorithms  using  a  linear  number  of 
processors  on  a  concurrent-read  concurrent-write  parallel  RAM. 


(a)  A  first  version  of  this  research  appeared  in  [Atallah]. 

(b)  This  research  was  supported  by  DOE  grant   DE-ACn2-76ER03077   and   by   NSF 
grant  NSF-MCS79-21 258 . 


/•/ 


I  .   I nt roduc  t  t on 

The  family  of  raodels  of  compvit^clon  used  In  this  paper  Is  the  prir.illel 
randora-access-machlnes  (PRAMs).  All  memhers  of  this  family  employ  p  synchronous 
processors  all  having  access  to  a  common  memory.  We  mention  three  members  of  the 
PRAM  family  In  descending  order  of  strength.  In  a  concur  rent- read 
concurrent-write  (CRCW)  PRAM,  simultaneous  reading  from  the  same  memory  location 
is  allowed  as  well  as  simultaneous  writing.  In  the  latter  case  the  lowest 
numbered  processor  succeeds.  A  concur  rent- read  exclusi ve-wr i te  (CRKW)  PRAM  allows 
simultaneous  reading  into  the  same  memory  location  but  not  simultaneous  writing. 
A.n  EREW  PRAM  does  not  allow  simultaneous  reading  or  writing.  See  |Vishkin83bl  for 
a  recent  survey  of  results  concerning  the  PRAM  family. 

In  the  present  paper  we  consider  the  problem  of  finding  Ruler  tours  in 
directed  and  undirected  Euler  graphs.  Parallel  algorithms  which  run  in  time  O(log 
|v|)  using  a  linear  number  of  processors  on  a  CRCW  PRAM  are  given. 

Consider  an  Euler  digraph  obtained  from  a  tree  by  replacing  each  one  of  its 
edges  by  two  anti-parallel  edges.  It  was  known  how  to  compute  an  Euler  circuit  of 
such  a  graph.  The  main  idea  of  Che  algorithm  is  quite  elegant:  We  find  a  spanning 
tree  of  some  auxiliary  graph  and  an  Euler  circuit  of  the  Euler  digraph  obtained 
from  the  tree.   Then  we  expand  this  Euler  circuit  to  an  Euler  tour  of  the  graph. 

II .   Basics 

The  following  definitions  and  propositions  follow  [Even].  Let  H(U,F)  be  a 
finite  undirected  graph.  An  Euler  path  of  H  is  a  path  e,,e2,...,eip|  such  that 
every  edge  appears  on  it  exactly  once.  An  undirected  graph  which  has  an  Ruler 
path  Is  an  Euler  graph . 

Proposl  t  ion  \^.         H  is  an  Euler  graph  iff  It  is  connected  and  exactly  two  vertices 
are  of  odd  degree  or  all  vertices  are  of  even  degree.   In  Che  latter   case,   everv 
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Kii  1  f  r  path  ot  Che  >^r;ii)li  Is  a  'Ircult  \i\<\     in  the  tonner  .asn  i-very  Kuler  path  leads 
from  one  of  the  two  odd  de)>rep  verloes  to  the  other. 

Let  r,(V,E)  be  a  finite  directed  graph.  An  F.u  I  e  r  path  of  C,  is  a  path 
e,  ,62  ,  .  .  .  ,  e  I  r.  I  such  that  every  edge  appears  on  It  exactly  once.  A  directed  graph 
which  has  an  Euler  path  Is  an  Euler  digraph.  The  underlying  undirected  graph  of  a 
digraph  is  the  graph  resulting  from  the  digraph  If  the  direction  of  the  edges  Is 
ignored.  We  use  deg  r(v)  (resp.  deg.  (v)  )  to  denote  the  In-degree  (resp. 
out-degree)  of  vertex  v. 

Proposition  2.  G  is  an  Euler  digraph  Iff  its  underlying  graph  is  connected  and 
one  of  the  two  following  conditions  holds: 

(1)  There  are  two  vertices  u  and  v  such  that,  ^^Squc^"^^  ~  deg^^^Cu)  +  1  and 
deg^^j.(v)  +  1  =  degj^j^(v),  while  for  all  other  vertices  w,  deg^^j^Cw)  =  deg|^j^(w). 
In  this  case  every  Euler  path  starts  at  u  and  ends  at  v. 

(2)  For  all  vertices  w,  deg  (^(w)  =  degj^pj(w).  In  this  case  every  Euler  path  is  a 
circuit . 

For  both  directed  and  undirected  graphs  we  restrict,  w.£.g.,  the  discussion 
to  Inputs  which  are  Euler  graphs  in  which  every  Euler  path  is  a  circuit.  Let  us 
justify  it:  (I)  Given  an  input  graph  we  can  check  within  the  efficiency  bound 
claimed  in  this  paper  if  Ic  is  an  Euler  graph,  if  it  has  Euler  circuits  and,  if  It 
has  Euler  paths  but  not  Euler  circuits,  what  are  the  endpoints  of  these  paths. 
For  this  apply  the  above  propositions.  The  only  non-trivial  detail  is  checking 
connectivity  of  underlying  undirected.  For  this  use  the  O(log  n)  time  algorithm 
of  [Shiloach  and  Vlshkin]  which  uses  a  linear  number  of  processors.  (2)  If  we  are 
given  an  Euler  graph  (or  digraph)  in  which  every  Euler  path  is  not  a  circuit,  It 
is  possible  to  "add"  a  single  edge  to  It  such  that  In  the  new  graph  every  Euler 
path  is  a  circuit.  After  finding  an  Euler  circuit  "delete"  the  added  edge  to 
obtain  an  Euler  path  in  the  original  graph. 
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Noce.  Wht'iiever  sfieclfylng  the  number  of  proci'ssors  vised  we  l/iuire  ihe 
constant  factor.  Namely,  Instead  of  "0(t)  time  using  0(p)  processors"  we  write 
"0(t)  time  using  p  processors".  Since  we  can  always  save  a  constant  factor  in  the 
number  of  processors  at  a  cost  of  the  same  constant  factor  in  the  rvinnlng  time. 

Subsection  III. I  gives  a  high-level  description  of  the  algorithm  fcir  finding 
Euler  circuits  in  digraphs.  A  detailed  description  follows  in  Subsection  III. 2. 
Section  IV  shows  how  to  extend  the  algorithm  of  Section  HI  for  finding  Ruler 
circuits  of  undirected  graphs. 

Ill .   Euler  Digraphs 

III .  1  A  high-level  description  of  Che  algorithm 

Step  1 .  Partition  Che  edges  of  che  InpuC  graph  G=(V,E)  inco  palrwise 
disjoint  circuits. 

Step  2_.  CorapuCe  an  auxiliary  undirected  bipartite  graph  G,=(V|,Ei)  defined 
as  follows.  G^  has  Cwo  seCs  of  verclces:  "clrculc-verc ices"  and  "real-vertices". 
Each  circuit  of  Scep  1  is  a  circuit-vertex  of  G..  Each  vertex  of  G  Is  a 
real-vertex  of  G,.  There  is  an  edge  beCween  a  real-vertex  and  a  circuit-vertex  If 
Che  (corresponding)  vercex  lies  on  Che  (corresponding)  circuit  in  G. 

Step  3.  Find  a  spanning  Cree  T=(V,,E,')  of  Gi.  Replace  each  edge  of  T  by  cwo 
anci-parallel  edges  (wlch  che  same  end-polncs)  co  obcaln  an  Euler  digraph  T. 

Scep  A_.  Find  an  Euler  clrculc  of  T  and  use  ic  Co  guide  che  sc  1  tching  of  the 
circuits  chac  were  found  In  SCep  1  Inco  an  Euler  circuit  of  G. 

Before  proceeding  to  a  detailed  descripclon  of  che  algorithm  we  summarize  a 
solution  of  [Tarjan  and  Vlshkin)  for  computing  an  Euler  circuit  of  T.  The  circuit 
is  computed  into  a  vector  FOLLOW.  For  each  vertex  v  of  V,  we  do  the  following. 
Let  deg(v)  =  d  in  T  and  let  the  d  adjacent  edges  of  v  in  T  be  ( v ,uq },...,( v ,u^_ ,) . 
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Asslgn  F()LL()W((u^,v))  :=  (v.u,^,  ,^^^,,  ^,  )  for  0  v  i    <      1-1. 

(This  Is  called  a  circular- success  ton  assignment ) 

FOLLOW  gives  for  each  edge  of  T  a  successor.   Thus,  given   an   edge   of   T   FOLLOW 

gives   a  path  that  starts  at  this  edge.   It  Is  easy  to  verify  that  FOLLOW  gives  an 

Euler  circuit  of  T.   Step  ■^.1  In  the  lower-level  description  that   follows   Is   of 

particular  interest.   There,  we  describe  how  the  circular  sucresslon  assignment  is 

actually  obtained. 

III. 2  The  algorithm 

Input  form.   The  vertices  of  G  are  represented  by  the  numbers   l,...,n.    The 
edges  are  represented  by  a  vector  of  length  |E|=m. 
We  employ  n  +  m  processors  In  the  CRCW  PRAM. 

Step  1 .  I .  Sort  all  edges  (l,j)  In  E  by  the  lexicographic  order  (defined  as, 
(a,b)  <  (c.d)  if  a  <  b  or  a  =  b  and  c  <  d).  Let  OUT-EDGE( v,k)  be  the  k-th 
outgoing  edge  of  v,  1  <  k  <  deg^^j.(v).  Sort  all  edges  (i,j)  in  E  in  the  reversed 
lexicographic  order  (defined  as  the  lexicographic  order  on  (j,l)  ).  Let 
IN-EDGE(v,k)  be  Che  k-th  incoming  edge  to  v,  1  <  k  <  deg.  (v). 

Implementation  Remark  1.  We  can  apply  the  sorting  algorithm  of  [Ajtai,  Komlos  and 
Szemeredi]  for  these  sorts.  This  algorithm  runs  In  O(log  ra)  time  using  m 
processors  and  0(m)  space  on  an  EREW  PRAM.  The  constants  In  front  of  log  m  are 
large.  It  is  possible  to  get  O(log  m)  time  (with  moderate  constants)  using  m 
processors  but  O(nin)  space  using  "orthogonal  trees"  (see  [Thompson]  for  more  (5n 
this  simple  algorithm).  We  avoid  elaborating  on  this.  A  third  posslbllty  is  to 
apply  the  randomized  sorting  algorithm  of  [ Re  1 f  and  Valiant).  It  runs  In  time 
Odog  m)  almost  surely,  using  m  processors. 

Step  1.2.   for  all  v  and  I  <  k  <  deg,^(v)  pardo 


snc(:(  in-ki)(;k(  v,k) )  :=  nirr-Ki)(;f-:(  v.k);  NKxr(  in-ki)(;k(  v,k))  :  =  out-i-:i)(;k(  v,k  ) 

Explanat Ion .  Step  1.2  completes  partitioning  the  edges  of  G  Into  circuit  as 
required  by  Step  1  above.  The  vector  SUCC  (and  the  auxiliary  vector  NEXT)  defines 
a  successor  for  each  edge  In  E  such  that  each  edge  In  E  has  exactly  one 
predecessor.  Therefore  the  paths  defined  by  SUCC  partition  the  edges  of  E  Into 
(disjoint)  circuits.  Our  goal  In  the  following  steps  Is  to  correct  SUCC  to 
contain  an  Euler  circuit.  At  this  stage  the  auxl 1 lary . vec tor  NEXT  Is  an  exact 
copy  of  SUCC.  The  reason  for  using  NEXT  Is  that  for  later  In  the  computation  we 
need  the  present  values  of  SUCC,  while  the  computation  of  Step  2.1  modifies  NEXT. 
Complexity .   0(1)  time,  ra  processors. 

Step  2.1.   Select  the  edge  of  minimum  lexicographic  order  of  each  circuit   to 
"represent"   it  in  the  auxiliary  graph  G,.   For  this  we  use  the  vector  D.  D((l,j)) 
Is  initialized  to  (i,j),  for  all  (i,j)  e  E.  The  iterations  below  use  Che  "doubling 
technique"  ( [Wyllie]  ). 
Apply  flog  ra|  iterations: 
for  every  edge  e  pardo 

D(e)  :=  MIN{D(e) ,D(NEXT(e))};  NEXT(e)  : =  NEXT(NEXT(e)  ) . 
D((i,j))   of   all   edges   of   the   same   circuit   have   now   the   edge   of  minimum 
lexicographic  order  of  the  circuit.   ( I  log  ml  iterations  sufficed   since   none   of 
the  circuits  has  more  than  ra  edges). 

An  edge  e  represents  the  circuit  it  lies  on  if  e=D(e). 
Complexity.   O(log  m)  time,  m  processors. 

Step  2.2.   Compute  the  auxiliary  bipartite  undirected  graph  G,(V],Ej). 
Vj  =  (v;  V  is  a  vertex  in  V  or  It  is  a  circuit  defined  by  SUCC( 
Recall   that   the  first  kind  of  vertices  are  called  real-vertices  while  the  second 


-6- 
are  called  circuit-vertices.  (Note .  In  our  writing;  wf  use  i\\o  fol  lowing; 
conventions.  We  do  not  always  distinguish  between  rea I -vert  Lees  ot  G,  and 
vertices  of  G.  We  also  do  not  distinguish  among  a  circuit  that  was  found  In  Step 
1.2,  the  corresponding  c Ircult- vertex  of  G,  and  the  edge  representing  this  circuit 
In  the  computation.) 

There  Is  an  edge  between  a  circuit-vertex  and  each  of  the  real-vertices  on  It: 
^1  -    i(v,w};  w  Is  a  circuit-vertex,  v  Is  a  real-vertex,  and  v  lies  on  the   circuit 
w  In  0} 

Let  {v,w}  be  an  edge  In  E.  where  v  Is  a  real-vertex  and  w  Is  a  circuit-vertex. 
Observe  that  the  vertex  v  may  appear  on  the  circuit  w  more  than  once.  At  each 
appearance  v  Is  entered  through  another  edge  of  the  form  (l,v)  (In  G).  Let 
CERTIFICATE({v,w} )  be  one  of  the  edges  of  the  form  (l,v)  on  circuit  w.  (The  reason 
for  the  name  is  that  it  certifies  the  edge  (v,w}). 

Implementation.  A  proper  application  of  sorting  (along  similar  lines  to 
Implementation  Remark  1)  would  readily  yield  the  graph  Gj  and  selection  of  edges 
for  CERTIFICATE({v,w}),  for  all  {v,w}  e  Ej.  Complexity.  The  graph  G^  has  0(n+Tn) 
vertices  and  O(n-hTi)  edges.  Hence  Step  2.1  takes  O(log  n)  time  using  n+m 
processors . 

Step  3.1.  A  spanning  tree  T  =  (V,,E[)  of  Gi  is  computed.  For  this  we  use  an 
adaptation  of  the  connectivity  algorithm  of  [Shiloach  and  Vishkin]  as  described  in 
[Tarjan  and  Vishkin].  Complexity.  O(log  n)  time  using  n+m  processors  on  a  CRCW 
PRAM.  ^ 

Step  3.2.  Obtain  the  digraph  T  from  T  by  replacing  each  edge  (v,w}  by  a  pair 
of  anti-parallel  edges  (v,w)  and  (w,v),  T  Is  an  Euler  digraph  since 
^^^In^''^  =  dsgout^^^  ^°"^  every  vertex  v  In  T. 
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Complex  1  ty .   Implementation  remark  1  boumis  Step   ) .  .'   i"   'fl'>K   n)   iliin'   ,ind   n 
processors . 

The  goal  of  Step  A. I  Is  to  compute  Into  the  vector  SUCC  a  large  circuit  which 
will  alternate  between  edges  of  G  and  edges  of  T  and  have  two  properties:  (1)  The 
edges  of  T  will  appear  on  It  In  an  order  of  an  F.uler  circuit  of  T.  (2)  The  edges 
of  G  will  appear  on  It  In  an  order  of  an  Euler  circuit  of  G.  We  concentrate  on 
achieving  the  first  property.  Interestingly,  the  second  property  will  follow 
without  further  care. 

Step  4.1.  Recall  the  Idea  of  finding  an  Euler  circuit  of  T  using  circular 
succession  assignments. 

It   Is   trivial   to   obtain   circular  succession  for  each  real-vertex  v  of  V,:  Let 
deg(v)  =  d  In  T  and  let  the  d  adjacent  edges  of  v   In   T   be   { v,uq} , . . . , ( v ,u j_i } . 
Set  SUCC((u^,v))  :=  (v.u^+i  ^^^  ^) ,    for  0  <  i   <    d-1. 
It   is  more  Interesting  to  obtain  circular  succession  for  each  circuit-vertex  w  of 

^1- 

Let  deg(w)  =  d  in  T  and  let  the  d  adjacent  edges  of  w  in  T  be  {  Vq  ,w}  , . . .  ,  { Vjj_^  ,w} . 

Recall   that   Vq,.,.,vj_^   are  real-vertices.   Let  (l^.v^)  be  CERTIFICATE( { v^ ,w} ) , 

for  0  <  a  <  d-1.   All  these  d  certificate  edges  appear  on  the  circuit  w  (In  G)   in 

some   circular   order.   The  key  idea  in  this  step  is  to  use  this  circular  order  in 

order  to  obtain  a  circular  succession  with  respect  to  T.   This  is  done  as   follows 

(Say  that  following  Step  1  SUCC(  (1^  ,v^ ) )  =  (v^,j^),  0  <;  a  <  d-1.   That  is,  ( v^  ,  j^ ) 

is  the  successor  of  (^^,v^)    in  the  circuit  w) : 

SUCC((v^,w))  :=  (v^.j^);  SUCC(  ( i„  ,  v^ ) )  :=  (w.v^).  for  0  <  c  <  d-1. 

We  just  set  the  successor  of  (v^,w)  to  be  ( v^ , j^ ) .   We  proceed  along  the  circuit  w 
(on   the   path  defined  by  SUCC)  till  we  hit  an  edge  (i^g.v^)  which  Is  a  certificate 
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edge  for  some  ori^e  { w ,  v„  }  i^  n  T.  (It  might  be  ( v^^  ,  )^  )  Itself).  SlICC  will  lead  us 
from  ('-g.vg)  to  (w,vg).  If  we  Ignore  the  edges  of  w  then  (w,vg)  Is  the  successor 
of  (v^,w).  Thus,  we  obtained  circular  succession  relative  to  c  I  rcul t- ver t Ices  . 
Therefore,  SUCC  contains  an  Euler  circuit  with  respect  to  T  as  required  (In 
property  (I)  above). 

Claim.   SUCC  contains  an  Euler  circuit  with  respect  to  G. 

Proof.  Let  (l,j)  be  an  edge  of  G.  we  show  that  It  Is  on  the  circuit  defined  by 
SUCC.  Following  Step  I  (l,j)  is  on  some  circuit  w.  Say  that  w  Is  a  circuit-vertex 
of  degree  d  in  T.  The  circuit  w  contains  d  certificate  edges  which  partition  the 
edges  of  w  Into  d  paths:  Each  such  path  starts  at  the  successor  of  a  certificate 
edge  and  ends  at  the  next  cetlflcate  edge  on  w.  (i,j)  must  lie  on  one  such  path. 
Each  cetificate  edge  must  Lie  on  the  circuit  defined  by  SUCC.  Hence,  (i,j)  must 
lie  on  the  circuit  defined  by  SUCC.  Also,  a  second  appearance  of  this  path  must 
follow  the  traversal  of  all  other  edges  of  G  or  otherwise  we  get  a  contradiction 
to  the  circuit  of  SUCC  being  an  Euler  circuit  with  respect  to  T.  • 
complexity.   0(1)  time,  n  processors. 

Finally,  we  need  to  "clean"  Che  circuit  of  SUCC  from  the  edges  of  T.  Observe  that 
there  are  no  more  than  two  successive  edges  of  T  in  the  circuit  defined  by  SUCC. 
Therefore,  Step  4.2  will  finish  the  Job. 

Step  4.2.   For  all  edges  e  of  G  perform  two  "doubling"  iterations: 

l_f_  SUCC(e)  is  an  edge  of  T 

then  SUCC(e)  : =  SUCC( SUCC(e) ) 
Complexity.   0(1)  time,  ra  processors. 

Complexity  of  the  whole  algorithm.   O(log  n)  using  n+m   processors.    Subject 

* 
to  Implementation  Remark  1  we  can  limit  the  space  Co  O(n-Hn)  as  well. 

>  >  I 
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IV.   Undirected  Graphs 

Let  H  be  an  Euler  undirected  graph.  We  show  below  how  to  orient  the  edges  of 
H  In  a  way  which  result  In  an  Euler  digraph.  This  solves  our  problem  since  now  we 
can  apply  the  algorithm  of  the  previous  section  to  this  digraph. 

Let  H'  be  a  directed  graph  obtained  from  H  by  replacing  each  edge  {u,v}  of  H 
by  two  anti-parallel  edges  (u,v)  and  (v,u).  Our  problem  Is  to  decide  which  one  of 
these  directed  edges  will  replace  {u,v}.  Let  v  be  a  vertex  of  H.  Say  that  Its 
degree  Is  d  where  d  Is  some  even  number.  Let  { v ,ui  }, (v ,U2 },••.,( v,uj }  be  the 
edges  adjacent  to  v.  We  partition  the  edges  of  H'  Into  circuits  as  follows. 
For  every  odd  number  1,  1  <  1  <  d,  assign 

NEXT((uj^,v))  :=  (v.u^^^);  NEXT(  (u^^^  ,  v) )  :=  (v,u^) 
Observe  that  the  circuits  defined  by  the  vector  MEXT  are  partitioned  into  pairs  as 
follows:  for  every  circuit  C  there  Is  another  circuit  C'  which  contains  exactly 
all  the  edges  which  are  anti-parallel  to  those  of  C.  Our  goal  is  to  delete  all 
edges  of  C  or  C  from  H'  for  every  pair  of  circuits.  For  each  e  e  H'  compute  into 
D(e)  the  minimum  edge  (i,j)  in  the  lexicographic  order  in  its  circuit.  (Apply  a 
computation  similar  to  Step  2.1).  Now  for  each  edge  (u,v)  in  H,  if  D((u,v))  < 
D((v,u))  then  (u,v)  is  selected  and  (v,u)  is  deleted.  For  each  pair  of  circuits 
this  results  in  selecting  all  edges  of  one  and  deleting  all  edges  of  the  other. 
We  indeed  get  an  Euler  digraph  since  the  fact  that  H  satisfied  Proposition  1 
implies  that  the  digraph  satisfies  Proposition  2.  This  computation  takes  O(log  n) 
time  using  m  processors  and  so  does  the  application  of  the  algorithm  of  the 
previous  section  to  the  Euler  digraph.   Again,  OCn+ro)  space  suffices. 

Conclusion 

Using  general  simulation  results  of  [Eckstein]  or  [Vlshkln83al  It  Is  possible 
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to  run  our  algorithms  on  an  ERF.W  PRAM  (or  the  stronger  CREW  PRAM)  In  time  OClog^n) 
using  a  linear  number  of  processors. 

The  only  place  where  our  algorithms  used  simultaneous  access   Co   the   memory 
location  Is  In  employing  the  Shlloach-Vlshkln  connectivity  algorithm  (Step  3.1). 

Recently,  [Awerbuch,  Israeli  and  Shlloach]  obtained  Independently  results 
which  are  similar  to  the  present  paper.  We  feel  that  our  solutions  for  both  the 
directed  and  undirected  case  are  simpler  than  theirs.  Theirs  and  ours  algorithms 
for  the  directed  case  are  pretty  similar.  However,  our  presentation  is  slightly 
simplified  due  to  the  observation  that  it  Is  possible  to  use  the  known  solution  of 
[Tarjan  and  Vishkin]  for  finding  an  Euler  circuit  in  a  tree  where  each  edge  is 
replaced  by  two  anti-parallel  edges.  We  use  it  both  for  the  intuitive  explanation 
of  the  algorithm  and  as  its  main  step.  In  the  undirected  case  they  claim  that  it 
is  not  that  easy  to  find  'Eulerian'  orientation  in  logarithmic  time.  However,  our 
solution  for  this  problem  Is  very  short  and  simple. 


-11- 

REFERENCES 


[Acallah]  M.J.  Atallah,  "Finding  Euler  tours  in  parallel",  Proc.  Seventeenth 
Annual  Conf .   on  Information  Sciences  and  Systems,  1983,  685-689. 

[Awerbuch,  Israeli  and  Shiloach]  B.  Awerbuch,  A.  Israeli  and  Y.  Shiloach,  "Finding 
Euler  circuits  in  logarithmic  parallel  time",  Proc.  Sixteenth  ACM  Symp.  on 
Theory  of  Computing,  1984,  249-257. 

[Ajtal,  Komlos  and  Szemeredi]  M.  Ajtai,  J.  Koral6s,  and  E.  Szemer^di ,  "An  0(n  log 
n)  sorting  network,"  Proc.  Fifteenth  ACM  Symp.  on  Theory  of  Computing 
(1983),  1-9. 

[Eckstein]  D.M.  Eckstein,  "Simultaneous  memory  access,"  Technical  Report  TR-79-6, 
Computer  Science  Department,  Iowa  State  University,  Ames,  Iowa,  1979. 

[Even]  S.  Even,  Graph  Algorithms,  Computer  Science  Press,  1979. 

[Reif  and  Valiant]  J.  Relf  and  L.J.  Valiant,  "A  logarithmic  time  sort  for  linear 
size  networks",  Proc.  Fifteenth  ACM  Symposium  on  Theory  of  computing,  1983, 
pp.   10-16. 

[Shiloach  and  Vishkin]  Y.  Shiloach  and  U.  Vishkln,  "An  O(log  n)  parallel 
connectivity  algorithm,"  J.   Algorithms  3  (1982),  57-63. 

[Tarjan  and  Vishkin]  R.E.  Tarjan  and  U.  Vishkin,  "An  efficient  parallel 
blconnectivity  algorithm",  Technical  Report  It69,  Computer  Science  Department, 
New  York  University,  New  York,  New  York,  1983. 

[Thompson]  CD.  Thompson,  "The  VLSI  complexity  of  sorting",  IEEE  Trans.  Comput.  , 
(December  1983). 

[Vlshkin83a]  U.  Vishkin,  "Implementation  of  simultaneous  memory  address  access  in 
models  that  forbit  it,"  J^  Algorithms  4  (1983),  45-50. 

(Vishkln83b]  U.  Vishkin,  "Synchronous  parallel  computation  -  a  survey".  Technical 
Report  //71,  Computer  Science  Department,  New  York  University,  New  York,  New 
York,  1983. 

[Wyllle]  J.C.  Wyllle,  "The  complexity  of  parallel  computation",  Technical  Report 
TR  79-387,  Department  of  Computer  Science,  Cornell  University,  Ithaca,  New 
York,  1979. 


v'c; 


This  book  may  be  kept 

FOURTEEN    DAY 

.  fi„e  .m  be  cba..ea  fo.  eacb  day  tbe  book  Ukep^^ve^ 


NYU/CoiTipSci/TR-13  4         c.l 

Atallah 
Finding  Euler   tours    in 
parallel. 


LIBRARY 

N.Y.U.  Courant  Institute  of 

Mathematical  Sciences 

251  Mercer  St. 
New  York,  N.  Y.    10012 


